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The master equation for the cluster-size distribution function is solved numerically to investigate 
the electron-hole droplet formation claimed to be discovered in the direct-gap CuCl excited by 
picosecond laser pulses [Nagai et ai, Phys. Rev. Lett. 86, 5795 (2001); J. Lumin. 100, 233 
(2002)]. Our result shows that for the excitation in the experiment, the average number of pairs 
per cluster (ANPC) is only around 5.2, much smaller than that (10 6 typically for Ge) of the well 
studied electron-hole droplet in indirect-gap semiconductors such as Ge and Si. 

PACS numbers: 71.35.Ee, 71.35.Lk 



In photoexcited semiconductors, electrons and holes 
are generated by the photon absorption and interact 
with each other via the Coulomb interaction. When the 
exchange- and correlation-induced attractions among the 
high density photogenerated carriers are strong enough to 
overcome the fermionic repulsive pressure and to reduce 
the chemical potential below the exciton energy, electron- 
hole droplets (EHDs) appear as first order transition. 1 
Investigation of the formation kinetics of EHD provides 
an opportunity to understand the quantum many-body 
effects. In semiconductor with an indirect gap and cor- 
respondingly long optical lifetimes, such as Ge and Si, 
EHDs have been investigated both theoretically and ex- 
perimentally in great detail. 2 Within the coexistence 
regime, relatively large EHD clusters are formed and the 
kinetics of the EHD formation is well understood in terms 
of the master equation 3,4 where the ANPC (n) is very 
large (for Ge (n) is about 10 6 ). 

Unlike indirect gap semiconductors, in direct gap semi- 
conductors the situation is very complicated because the 
system cannot reach quasiequilibrium during the carrier 
lifetime (typically in the order of 1 ns). Despite exten- 
sive experimental efforts, no decisive evidence of the EHD 
formation has yet been reported for a very long time. 5 
EHD formation has been predicted theoretically in po- 
lar CdS and CdSe. 6 Also, EHD-like behavior has been 
observed in them experimentally in which time-resolved 
emission measurements are performed. 7 ' 8 This method 
has also been applied to GaAs. It reveals that electron- 
hole plasma created by band-to-band excitation has ex- 
tremely high temperature due to the overheating caused 
by the reduction of band gap. Very recently, Nagai et 
al. used CuCl which has large binding energy and band 
gap (binding energy of IS-Z3 exciton E e2; =213 meV and 
band-gap energy E g =3.396 eV) to overcome the problem. 
They therefore are able to study the formation of EHD in 
CuCl via time-resolved emission and pump-probe reflec- 
tion measurements by femtosecond excitation in the mid- 
infrared region under resonant exciton excitation. They 
reported that clear EHD-like behavior in CuCl has been 




FIG. 1: Cluster concentration versus time and number of e- 
h pairs per cluster for CuCl at T = 4.2 K under Gaussian 



pulse excitation, t r 
po = 2.0 x 10 20 cm 
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However, due to the fact of the short lifetime of the car- 
rier, one expects only small clusters can be formed and 
sustained. Two decades ago, Haug and Abraham per- 
formed a theoretical investigation of the EHD formation 
in GaAs. 12 They reported for the typical picosecond exci- 
tation the ANPC is at most around 2 where there should 
be no EHD formation can be found. In this paper we 
use the master equation theory developed by Haug and 
Abraham 12 to investigate the femtosecond excitation in 
CuCl in the experiments of Nagai et ai. 10,11 This the- 
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ory has been proven to be very sucessful in dealing with 
the EHD formation in indirect gap semiconductors. 3,4 In 
this study, we focus on how big the ANPC can be in this 
material and how this number is affected by different pa- 
rameters. 
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FIG. 2: Concentrations of selected clusters versus time for 
CuCl at T=4.2 K under Gaussian pulse excitation as in figure 
1. 

The concentration f n (t) to find at a given time t a 
cluster with n e-h pairs per unit volume is given by the 
following master equation equations: 12 For n > 1, 

dfn . m 

q_I_ ~ Jn— 1 Jn {i-/ 
]n = 9nfn - ln+lfn+1 (2) 



e„ is given by the Richardson-Dushman current 3 

f m*kT\ 3/2 g n _i -i/3\m 
En ^ 1 ) —jr-cMi-fi + c™ ' )® (7) 

7 is the degeneracy of the exciton level, <fr is the binding 
energy of the liquid with respect to the exciton. ccm -1 / 3 
is the correction of the binding energy with a denoting 
the surface energy of the liquid. Typically we select 7=1. 
a = Cipporo with C ~ 1, <p standing for the condensate 
energy per electron-hole pair and r = (S/iirpo) 1 ^ 3 . 1 The 
change of exciton (corresponding to n = 1) is treated 
separately: 

**=G(t)-5:f -iyi-5> (8) 

n=l n n=2 

where G(t) is the generation rate of excitons. 
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where j n is the net probability current between the clus- 
ters with n and n+1 c-h pairs, g n is the gain rate of a 
cluster with n c-h pairs and is approximated by 

g n = inRlvofi (3) 
v = ^kT/2nm* (4) 

« = Y RnPo (5) 

R n is the radius of the cluster, w is the thermal exciton 
velocity, po is the liquid density. The loss rate l n is the 
sum of the evaporation rate e n and the recombination 
rate n/r n : 

l n = e n + n/T n . (6) 



FIG. 3: ANPC versus time for CuCl at T=4.2 K under Gaus- 
sian pulse excitation. t r =0.2 ps, Go = 6.1 x 10 33 cm _3 s _1 , 
po = 2.0 x 10 20 cm~ 3 . Solid curve: r = 1 ns; Dashed curve: 
r = 0.5 ns; Dot-dashed curve: t—2 ns. 

We solve above equations numerically for the exam- 
ple of CuCl. The parameters for CuCl at 4.2 K are 
TO*=2.0m and (f> = 5.12 x 1CT 14 erg. 13 The value of 
the binding energy <j> is taken to be the biexciton bind- 
ing energy which is the lower limit of the binding energy. 
Our computation indicates that the result is almost un- 
changed if one increases <j) by one order of magnitude. 
We obtain p from r s = ^S/Annal, where a is the Bohr 
radius and r s =1.8. 10 a — 0.54 erg-cm -2 . In our calcula- 
tion we assume t„ = r. According to the experiment, we 
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choose a Gaussian pulse: 

G(t) = G oe - (t -* o)2/t ' 



(9) 



where the pulse width t r = 0.2 ps, and the peak genera- 
tion rate is Go = 6.1 x 10 33 cm _3 s _1 , corresponding to a 
excitation density of 4.3 mJ-cm~ 2 . The main results of 
our calculation are plotted in Figs. 1-3. 

Figure 1 shows the development of the distribution of 
cluster distribution. For the short pulse, the exciton con- 
centration fi(t) rapidly reaches its maximum near the 
pulse center, to = 2 ps. Then it decreases quickly due 
to the formation of the clusters. It is also noted from 
the figure that within 1 ps after the center of the pulse, 
clusters of all sizes reaches to their maximum values. As 
any cluster with size n > 2 comes from the exciton, f\ (t) 
sharply drops 5 orders of magnitude from its peak within 
that period of time. After that, because of the recombi- 
nation and the evaporation, almost all clusters start to 
decay, except for some small clusters (say n = 1 and 2) 
which increase a little bit due to the evaporation from 
bigger clusters. The lifetime of large cluster (say n = 50) 
is much shorter than that of the small one as the later 
gets compensation from the decay of the larger clusters 
and also the recombination rate n/r n is propotional to 
n. After 10 ns almost all cluster concentrations are below 
10 10 cm~ 3 . This can be seen more quantitatively in Fig. 

2 where the concentration of some selected clusters are 
plotted as functions of time. 

In order to see in average how large the cluster can be 
under the experiments by Nagai et al., we plot in Fig. 

3 the ANPC (n) as function of time for three different 
lifetimes r. (n) can be calculated from 
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J2n=l H fn 
X)n=l fn 



One notices from the figure that for the excitation in 
the experiment, ANPC is only 5.2, orders of magnitude 
smaller than the number in the EHD of indirect semicon- 
ductor such as Ge and Si. This result is consistent with 
the regular expectation of the excitation in the direct gap 
semiconductors. Moreover, if one assumes that the larger 
clusters are stable and therefore r„ (n > 2) is larger than 
T\. We demonstrate in the same figure the result when 
Ti = 1 ns but t„ = 10 ns when n > 2. One finds that this 
does not affect the ANPC, although the metastable time 
is much longer. We further find that if we increase the 
excitation Go by one order of magnitude, we can only 
double the ANPC. These results indicate that the fea- 
ture of low ANPC is very robust. These results indicate 
that what measured in CuCl by Nagai et al. 10,11 may not 
come from the EHD formed from exciton gas, but instead 
possibly come from some bubbles of excitons in metallic 
liquid. 11 
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